Can a frustrated spin-cluster model describe the low-temperature physics of NaV205 ? 
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Recent experimental evidence suggest the existence of 
three distinct V-valence states (V^"*, V^^'^ and V"*"^) in the 
low-temperature phase of NaV2 05 in apparent discrepancy 
with the observed spin-gap. We investigate a novel spin clus- 
ter model, consisting of weakly coupled, frustrated four-spin 
clusters aligned along the crystallographic b-axis that was re- 
cently proposed to reconcile these experimental observations. 
We have studied the phase diagram and the magnon disper- 
sion relation of this model using DMRG, exact diagonalization 
and a novel cluster-operator theory. We find a spin-gap for 
all parameter values and two distinct phases, a cluster phase 
and a Haldane phase. We evaluate the size of the gap and 
the magnon dispersion and find no parameter regime which 
would reproduce the experimental results. We conclude that 
this model is inappropriate for the low-temperature regime of 
NaVaOs. 

PACS numbers; 75.10.d, 72.15.N, 71.20.B, 75. 10. J 
Introduction Recent investigations of electronically 
quasi one-dimensional (ID) transition metal compounds 
probe the limits of our understanding of the interplay 
between structural and electronic effects in such low di- 
mensional materials. InNaV2 05, a prototypical example 
for this class of materials, V-ions are arranged in ladders 
along the crystallographic b-direction. Measurements of 
the magnetic susceptibility in the high-temperature 
phase indicate the presence of only one equivalent V-sitc 
with valence V"'"'''^, consistent with a model where 
the electrons in bonding V-O-V orbitals along the rungs 
of the ladder form a ID Heisenberg chain . 

At Tc = 34 K the unit-cell doubles along the a- and b- 
and quadruples along the c-axis in a phase transition 
of as-of-yet unknown origin. At the same time a spin- 
gap of Ami„ = lOmeV opens 0] and charge ordering 
2Y+4.5 ^ y+4 ^ y+5 gg^g rpj^g observed charge 

ordering is inconsistent with a generic spin-Peierls sce- 
nario Ql and raises the question about the driving force 
(lattice, magnetic or Coulomb) responsible for this tran- 
sition. Since NaV2 05 is an insulator, the discussion of 
the material is simplified by the introduction of pseu- 
dospins for the charge degrees of freedom that couple to 
the spin-degrees of freedom [0-0. The effective spm- 
Hamiltonian depends, consequently, on the pattern of 
charge order 111]! and may differ in the high- and the 
low-temperature phase. 

The occurrence of two well defined magnon-branches 
for T < Tc in NaV205 along the a direction (perpen- 



dicular to the chains), as measured by neutron scat- 
tering had been explained tentatively by a model, 
where the charge orders in a 'zig-zag' pattern in the low- 
temperature phase |]ri|. This proposal has been ques- 
tioned by recent analysis of the low-temperature crystal 
structure [|l2yi3| ]. Based on bond-charge models, the exis- 
tence of three different V-valence states (V"*"*, V+'*'^ and 
V+^) has been proposed ||l^,|l^, as illustrated in Fig. 
|l|. In this analysis, pairs of V+^'^ form dimerized spin- 
chains on every other ladder, which alone could explain 
the observed spin-gap A puzzle is posed however, 

by the presence of free isolated moments on the V"'"'' ions 
on the remaining ladders, which is inconsistent with the 
existence of such a gap. 

As one possible reconciliation, Boer e< al. ||l^ recently 
proposed that the V^^ moments are quenched by their 
interaction with the neighboring V+^'^ sites of the adja- 
cent dimerized V-O-V ladder. Within this model, clus- 
ters of six Vanadiums each (and with four spins) would 
be weakly coupled and the observed spin-gap would arise 
not from the dimerization but locally from the gap of the 
isolated clusters. 

To distinguish between these fundamentally differ- 
ent mechanisms we study this model by a series of 
complementary approaches, using DMRG exact- 
diagonalization and a novel bond-cluster theory, to map 
all physically relevant regions of its phase diagram. We 
find that the ground state varies continuously from a 
cluster-phase for large J' to a Haldane-phase for small 
J' (see Fig. |l|). We evaluate the gap and the dispersion 
and find that there is no parameter regime that would 
explain the neutron-scattering data |l^ ]. 

The spin-cluster modelWe denote by S„^i (i = 1, . . . , 4) 
the four spins of the n'th cluster, compare Fig. I. The 
Hamiltonian is then 

H — Ji S„_i • S„_2 + J2 S„_i • S„+i_2 



+ J'Y. + Sn,2) • (S„^3 + S„,4) , (1) 

n 

where Ji J(l -I- 5) and J2 = J(l - 5) (with Ji, J2, J' > 
0). S is the degree of dimerization. For J' — the S„ 1/2 
form a dimerized chain with an in-chain gap ~ J6^/^ . A 
particular property of Eq. (||) is the local coupling to the 
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total spin S„_3 + S„^4, which is consequently a (locally) 
conserved quantity, {Sn,3 + S„^4)^ = S'„(S'„ + 1) for any 
n. In the ground-state S'„ = 1. 



• 0(1) 




FIG. 1. The spin-cluster model for NaV205. On the left 
the charge valency in one V-0 plane obtained from two crystal 
structure determinations On the right the proposed 

cluster spin-model Note, that two V^*'^ on one rung 

share one electron. 

We consider first an isolated cluster and denote by Sy 
and tfj the wavefunctions of the singlet and of the triplets 
{a = —1,0, +1) of the spins i and j = 1, . . .4,). The 
low-energy states are 

V"! = [^12^34 ~ ''12^34 ^ ''12^34] ' (2) 

tp2 = S12S34, ■03 = Sl2*34 > (3) 

^4° - ^ [te - ^r244] , (4) 

where ip^ is the — component of the triplet ip2 ■ The 
corresponding energies are Ei = —2J' + Ji/ A, E2 ~ = 
-3Ji/4 and E^ = -J' + Ji/4 

For J'/Ji > 0.5 the singlet 01 is the ground-state (we 
denote this region the 'cluster phase'). For J'/Ji < 0.5 
the ground-state of the isolated cluster is four-fold degen- 
erate, the singlet ip2 and the triplet ip^ have the same en- 
ergy. Note that the intercluster coupling J2 will not mix 
ip2 and ips, since the local spin S„^3 + S„^4 is conserved. 
Intercluster coupling J2 will lead to an antiferromagnetic 
interaction Jh ~ {J' J2)^/Ji between the moments of 
the ■03 states, as can be evaluated easily in second-order 
perturbation in J2 (using the complete set of eigenstates 
of the cluster). The total energy is therefore lowered by 
J2 when all cluster states are 03. The 5* = 1 moments 
of the ips states thus form an effective spin-1 chain with 
a Haldane gap = 0.41050 J// Q. We denote this 
region therefore the 'Haldane phase'. 



We have evaluated the energy gap of the spin-cluster 
model by DMRG |lj] , using the finite-size algorithm with 
open boundary conditions for systems with L = 32 and 
L = 64 spins. The ground-state has iVf = L/2 up-spins 
and Ni — L/2 down-spins. We retained typically 60 
states of the density matrix, checking the convergence by 
additional calculations with 40 and 90 states respectively. 
We evaluated the gap by two complementary methods, 
namely (i) by targeting two states in the sector with 
iV| — L/2 = Ni and (ii) by targeting the ground states in 
(a) the sector with = L/2 = Ny and (b) = L/2 + 1 
and Ni — L/2 — 1. We find complete consistency and 
present the results in Fig. || for some selected values for 
the dimerization 5. The finite-size corrections are smaller 
than the symbol sizes. We find a rapidly decreasing gap 
as a function of decreasing J'/ J and a smooth crossover 
between the cluster- and the Haldane phase. As the sym- 
metry of these two phases is the same, we do not expect 
a phase transition in the thermodynamic limit. 

J' /J, 




J'/ J 

FIG. 2. The (singlet-triplet) gap in units of J for fixed 
(5 = and (5 = 1/3 as a function of J'. Filled Symbols: 
DMRG data. Lines: Results from the cluster operator theory. 
Inset: The three low- lying energy levels Ei, Ez and E4, for an 
isolated cluster in units of Ji as a function of J' / Ji- The 
corresponding wavefunctions ip^ and ip'i are triplets, ipi is a 
singlet. 

Cluster- operator theory In the cluster-phase two low- 
lying triplet modes, ■03 and ^04, are relevant. In order to 
take the effect of the intercluster coupling J2 into account 
we describe the seven degrees of freedom of cluster n by 
bosonic degrees of freedom: for the singlet (-01) fojj 3 

and b]^^^ for the triplets (1/^3 and ^04). The low-lying 
singlet i/'2 does not couple and may be disregarded here. 
This approach generalizes the bond-operator theory for 
dimcrized spin-chains Eol to the case of spin-clusters. 
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The constraint s+s„ + Er,a &l,r,a^n,T,a = 1 (r = 3,4) 
restricts the bosonic Hilbert space to the physical one. 
The spin-operators take the form 



'-'n.l 



/2 



± 



12 



V6 



(5) 



Note, that there are no terms ~ b'l^,T,cfin r' a' 
responding to triplet-triplet interactions. In linearized 
Holstein-Primakov approximation (LHP), we substitute 
sjj ^ 1 and s„ ^ 1 in Eq. (||) and in similar expressions 
This approximation retains spin-rotational in- 
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variance and we may disregard the index a = — 1, 0, 1 for 
the triplet operators. We obtain for the LHP Hamilton- 
operator in momentum space 

with Ho = J2k,r ^-^I,r^,r i^r = - Ei) . Thc intcr- 
cluster coupling is given by 
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analytic 
L = 16, exact 
L = 64, DMRG 



J' = 2.0, 8 = 0.2 
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FIG. 3. The magnon dispersion for J' = 2J,5 = 0.2. 
Note the zero of energy. The lines are the result of the clus- 
ter-operator theory, the circles of an exact diagonalization 
study with 16 spins (periodic boundary conditions). The cross 
denotes the DMRG result. Inset: The value of J as a func- 
tion of dimerization 5 needed to fit the measured dispersion 
of NaV205 (LHP result, for J/2 < J' < J). 



-z2V2 sin(25fc) (foi 4 - 6^6, 3 



(6) 



and 



-i2V2 sin(2fefc) ft^^^^.g + h.c. 



(7) 



— analytic 
• L = 16, exact 
X L = 64, DMRG 




J' = 0.8, 8 = 0.2 
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FIG. 4. The magnon dispersion for J' = 0.8 J,S = 0.2.The 
symbols as in Fig. H. Note that the cluster-operator theory 
overestimates the magnon-dispersion. 

Here b — 3.611 A is the lattice constant of the high- 
temperature phase. Note the opposite sign in the dis- 
persion of two triplets. It is straightforward to diag- 
onalizc H^^"P\ We define c = ( J2/6) cos(2&fc), 2t = 
A| + A§ + 2c(2A4 - A3) and s = A| A| + 2cA3A4(2A3 - 
A4) — 2A3A4 j|/9. The dispersion = iiJ±{k) of the two 
magnon branches (each branch is three-fold degenerate) 
in LHP-approximation is then 



col 



(8) 



We have included the results for the magnon gap in 
Fig. ||. For large ratios J'/ J the LHP-result becomes 
asymptotically exact, in this limit it is equivalent to per- 
turbation theory in J2. In the LHP-approximation the 
transition to the Haldane phase is signaled by a vanish- 
ing of the energy gap, the crossover cannot be described 
by the cluster-operator theory. 

In Fig. ^ we present the magnon-dispersion Eq. (|^) 
for J' = 2J and compare the LHP-results (lines) with 
an exact-diagonalization study of a system with 16 sites 
(filled circles) |2^. The agreement is very good, due 
to the large gap and (correspondingly) small correlation 
length. Note that the low-lying magnon, which corre- 
sponds to V'4 (see inset of Fig. H), has its minimum at 
k = 7r/(2&). 

In Fig. ^ we present the magnon-dispersion Eq. (|^) for 
J' — 0.8 J which is closer to the transition to the Haldane 
phase. The agreement with the exact diagonalization and 
the DMRG data is not good, since the precursors to the 
Haldane phase are not included in the cluster-operator 
theory. The low-lying magnon, which corresponds to 
(see inset of Fig. ^), has its minimum now at fc = and 
k = n/b and a maximum at fc = 7r/(26), as measured 
by neutron scattering p6| . The cluster-operator theory 
substantially overestimates the size of the magnon dis- 
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persion relative to the exact-diagonalization result near 
to the Haldane phase. The physical reason for this dis- 
crepancy can be understood: The lattice constant of the 
effective spin-1 chain in the Haldane phase is 2b and the 
minimum of the magnon dispersion is therefore at tt/ (26) 
in the Haldane phase |2^. It changes therefore at the 
crossover from the cluster phase and the Haldane phase. 
This change in the location of the gap is not included in 
the cluster-operator theory. 

Discussion The exchange constant along 6 is J ~ 
529 — 560 K P,p2| in the high-temperature phase of 
NaV205 and the inter-ladder coupling is probably very 
small, a J'/J ~ 1/45 has been found in an analysis of the 
magnon-dispersion for T < Tc in a model with zig-zag 
charge order This small ratio is consistent with the 
very small coupling along a found in a LDA-study [Q. 
There are, however, two reasons why J' might be larger 
in the low-temperature phase, (a) As noted by Mack and 
Horsch , there is a near cancellation for T > Tq in be- 
tween paths with intermediate singlet and triplet states 
and energies Eg ft: J' — t1y{l/Es — l/Et), where t^y is 
the V-V hopping matrix element in a direction. A cor- 
responding calculation for T < Tc in the phase shown 
in Fig. |l| yields J' = 2tly{l/Es + 1/U - l/Et) {U is 
the onsite Hubbard-JJ). (b) txy might be substantially 
larger in the low-temperature phase, since the smallness 
of tncy for T > is a subtle band-structure effect [Q . We 
have therefore scanned the complete phase diagram of the 
spin-cluster Hamiltonian in order to determine whether 
there exists a parameter range able to fit the neutron 
scattering data. 

We have tried to reproduce, within the spin-cluster 
model, four known properties of NaV2 05: (i) The gap 
(averaged over ka) is Amin = lOmeV. (ii) The max- 
imum of the dispersion of the lowest magnon branch 
is at 7r/(26), the minimum at and Ti/h. (iii) The 
value of the maximum of the dispersion of the lowest 
magnon branch is ls.max ~ 40meV [^,|3|, i.e. the ratio 
is Amaa;/A„ij„ « 4. (iv) The value of the coupling along 
6 is J « 441 K = 38 meV for T < Tsp 

Condition (ii) implies that only the cluster-phase of 
Hamiltonian Eq. (^) with J' < Ji is a candidate for 
the low-temperature phase of NaV205. This implies 
Ji/2 < J' < Ji. Within the cluster-operator theory 
one obtains Amax/^min = 4 for values of J' near to 
the gap closing. One needs consequently large coupling 
constants J (see inset of Fig. in order to reproduce 
Amin = lOmeV. We have evaluated the values of J' 
and J needed to reproduce the gap-ratio as a function of 
dimerization S and find a minimum in J for 6 — 0.2 (see 
inset of Fig. |). This minimum is J « 126meV, substan- 
tially larger than the experimental value J « 38 meV. 
Note, that the cluster-operator theory overestimates the 
dispersion in this phase and underestimates the value 
of J needed. We therefore conclude safely, that the model 
is not able to reproduce the measured magnon-dispersion 



of NaV205 and that Eq. (|i|) is unlikely to be the appro- 
priate model for the low-temperature phase of NaV205, 
at least in its one-dimensional version. It might be possi- 
ble, in principle, that two-dimensional couplings change 
the scenario obtained in the present study, though we 
note, that an increase in dimensionality does, in general, 
reduce the size of a spin-gap. 

We would like to acknowledge discussions with P. Lem- 
mens and the support of the German Science Foundation, 
the BMBF and the Fonds der chcmischcn Industrie. 



[1 
[2: 

[3: 
[4: 

[5 
[G 

[7 

[8 

[9: 

[10: 
[11 

[12 
[13 
[14 

[15 

[16 
[17 



[is: 

[19 



[20 

[21 
[22; 

[23' 
[24; 



H. Smolinski et al, Phys. Rev. Lett. 80, 5164 (1998). 
H.G. von Schnering et al, Z. Kristallg. 213, 246 (1998); 

A. Meetsma et al., Acta Cryst. C 54, 1558 (1998). 

P. Horsch and F. Mack, Euro. Phys. J. B 5, 367 (1998). 
M. Isobe and Y. Udea, J. Phys. Soc. Jap. 65, 1178 (1996). 
Y. Fujii et al, J. Phys. Soc. Jap. 66, 326 (1997). 
T. Ohama, H. Yasuoka, M. Isobe and Y. Udea, Phys. 
Rev. B 59, 3299 (1999). 

H. Seo and H. Fukuyama, J. Phys. Soc. Jpn. 67, 2602 
(1998). 

P. Thalmeier and P. Fulde, Europhys. Lett. 44, 242 

(1998) . 

M. Mostovoy and D. Khomskii, Sol. St. Comm. 113, 159 

(1999) . 

Debanand Sa and C. Gros, preprint. 

C. Gros and R. Valentf, Phys. Rev. Lett. 82, 976 (1999). 

J. Liidecke, A. Jobst, S. van Smaalen, E. Morre, C. Geibel 

and H.-G. Krane, Phys. Rev. Lett. 82, 3633 (1999). 

J.L.de Boer, A.M. Meetsma, J. Baas and T.T.M. Palstra, 

preprint. 

S.R. White, Phys. Rev. Lett. 69, 2863 (1992). ; L Peschel 
et al. (eds). Density -Matrix Renormalization (Springer, 
Berlin, 1999) 

S. van Smaalen and J. Liidecke, Europhys. Lett. 49, 250 
(1999). 

T. Yosihama et al, J. Phys. Soc. Jap. 67, 744 (1998). 
The remaining eigenstates are the S — 2 state -05 
(/3 = — 2, . . . , 2) and the triplet -i/jg = ^12^34, with en- 
ergies E5 ^ J' + Ji/4 and Eq = Ji/4. 
S.R. White and D.A. Huse, Phys. Rev. B 48, 3844 (1993). 
S. Sachdev and R.N. Bhatt, Phys. Rev. B 41, 9323 (199); 

B. Normand et al., Phys. Rev. B 56, R5736 (1997); W. 
Brenig, Phys. Rev. B 56, 14 441 (1997). 

For the method used see V.N. Muthukumar et at. Phys. 
Rev. B 54, R9635 (1996). 

See, e.g., I. Affleck, J. Phys. Cond. Matt. 1, 3047 (1989). 
F. Mila, P. Millet and J. Bonvoisin, Phys. Rev. B 54, 11 
925 (1996). 

T. Yosihama, private communication. 

M. Weiden et al, Z. Phys. B 103, 1 (1997). 



4 



